Quark and Lepton Masses from Gaussian Landscapes 
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The flavor structure of the Standard Model might arise from random selection on a landscape. 
We propose a class of simple models, "Gaussian landscapes," where Yukawa couplings derive from 
overlap integrals of Gaussian wavefunctions on extra-dimensions. Statistics of vacua are generated 
by scanning the peak positions of these zero-modes, giving probability distributions for all flavor 
observables. Gaussian landscapes can broadly account for all observed flavor patterns with very 
few free parameters. For example, the generation structure in the quark sector follows from the 
overlap integrals for both the up and down type Yukawas sharing the localized wavefunctions of 
the quark doublets and the Higgs boson. Although Gaussian landscapes predict broad probability 
distributions, the flavor observables are correlated and we show that accounting for measured flavor 
parameters creates sharper distributions for future neutrino measurements. 



1. Introduction: The Standard Model (SM) has en- 
joyed remarkable success at explaining laboratory data. 
Nevertheless, it requires 28 parameters to be set by hand, 
thus begging for a more fundamental description. Most 
of these parameters appear in the flavor interactions 
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where q, I (u, d, e, V) are the left (right) handed quark 
and lepton fields, h is the Higgs boson, and <j> represents 
the additional scalar(s) responsible for giving Majorana 
masses to the right-handed neutrinos v. Motivated by the 
success of unified gauge symmetries in describing the SM 
gauge couplings, the conventional wisdom is that some 
flavor symmetry is behind the patterns seen in the 22 
flavor observables stemming from Eq. JT]). Indeed, it is 
well known that mass hierarchies and small mixing angles 
can arise from small flavor-symmetry breaking parame- 
ters when different SM generations feel different levels of 
the flavor-symmetry breaking On the other hand, 
we have yet to discover any precise, compelling relations 
among the flavor observables that would confirm a fun- 
damental symmetry principle. 

The cosmological dark energy may be evidence for a 
huge landscape of vacua, with the observed value of the 
cosmological constant resulting from environmental se- 
lection for large scale structure Q. The current under- 
standing of string theory seems to support the existence 
of such a huge landscape of vacua. For example, in the 
Type IIB string description flux compactification gener- 
ates a statistical distribution of complex structure mod- 
uli [H, which then determines the statistics of Yukawa 
couplings. Thus, some kind of statistical randomness 
may be involved in the Yukawa couplings and this ran- 
domness may explain the absence of relations among 
masses, mixing angles, and CP phases that would other- 
wise reflect some fundamental symmetry principle. This 
philosophy has been pursued in Refs. [J, Q. 

Can a theory of flavor using sheer randomness explain 



the various qualitative patterns among the flavor observ- 
ables? We consider seven major features of flavor to be i) 
the hierarchical masses in the quark and charged lepton 
sectors, ii) the pairing structure (i.e. small mixing an- 
gles) of the quark sector, iii) the generation structure of 
the quark sector (that is, the electroweak pairing between 
the two heaviest, middle and lightest quarks), iv) the ab- 
sence of pairing structure (i.e. large mixing angles) in the 
lepton sector, v) the mixing angle #13 of the lepton sector 
being not so large as the other lepton mixing angles, vi) 
the hierarchy among Yukawa eigenvalues being largest in 
the up-quark sector, and vii) the CP phase in the quark 
sector being of order unity. Assuming simple statistical 
behavior of Yukawa couplings, Ref. [5[ describes features 
(i) and (ii), while Ref. [4[ describes feature (iv). 

Although these results are encouraging, they have sev- 
eral shortcomings. What is often broadly referred to as 
the generation structure of the SM includes features (i)- 
(iii), and is more confounding in light of the large- mixing 
neutrino oscillations. Thus we consider it crucial to ob- 
tain features (i)-(iv) all at once in a single, simple frame- 
work. Refs. [i[ and @ use different schemes to describe 
the charged fermion and neutrino sectors, and fail to ac- 
count for feature (iii). Furthermore, these models lack a 
compelling motivation in terms of a more fundamental 
theory. We introduce "Gaussian landscapes" as models 
of subsets of the landscape expected from compactifica- 
tion of Heterotic string theory. As in Ref. [|| , mass hier- 
archies, feature (i), arise without flavor symmetries due 
to small overlap integrals involving Gaussian wavefunc- 
tions on extra dimensions. Whereas in Ref. Q there is 
no landscape and these wavefunctions are positioned by 
hand to agree with data, we take the peak positions to 
scan randomly over the geometry of extra dimensions, 
with each possibility corresponding to a distinct vacuum 
of the landscape. Treating our universe as a typical, ran- 
dom selection from such a landscape, we find that Gaus- 
sian landscapes can broadly account for each of the seven 
features of flavor described above. Ref. Q gives a more 
extensive account of the material in this letter. 
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2. Gaussian Landscapes: In the compactification of 
the Heterotic string theory, Yukawa couplings are calcu- 
lated by overlap integration of zero-mode wavefunctions 
over a six-dimensional compactification manifold. When 
the overlap of wavefunctions happens to be small, the 
relevant Yukawa couplings are small. It is known that 
the wavefunctions are approximately Gaussian on a base 
manifold B when a six-dimensional manifold is a torus 
fibration and the size of torus is small relative to B Q. 
Meanwhile, if the gauge field moduli vary from one part 
of the universe to another, then the corresponding zero- 
mode wavefunctions vary, and likewise the masses and 
mixing angles. The statistics of gauge field moduli, 1 
therefore, generate the statistics of flavor observables. 
Statistics generated in this way are basis independent. 

Aiming to extract the essence of string theory compact- 
ification with respect to flavor, we introduce the Gaus- 
sian landscape as a simplified version of the landscape. 
The Gaussian landscape posits that low-energy degrees of 
freedom have localized (Gaussian) zero-mode wavefunc- 
tions in some geometry of extra dimensions, which is an 
analogue for the base manifold B. In this letter we first 
use S 1 as the geometry, and briefly discuss the impact of 
the choice of geometry later. The wavefunctions 

tfiy; yf) ~ A cxp [-(1 + ir)(y - y?) 2 /(2<£)] (2) 

are taken to be complex with a phase related to the uni- 
versal parameter r. Here a = q,u,d,l,e,V,h,4> labels 
particle species and i the generation of fcrmions. The 
factor Ai is a normalization factor chosen so that 

M 5 [ dy\<p?(y;y?)\ 2 = l, (3) 
Jo 

where M5 is the cut-off scale of the four-dimensional ef- 
fective theory and L is the circumference of S 1 . In fact, 
the wavefunction (|2|) is made periodic on S , while main- 
taining the normalization in . Yet as long as the width 
d a is paramctrically smaller than the circumference L, 
the wavefunction is almost Gaussian. 

In the Gaussian landscape, the up-type Yukawa matrix 
derives from the overlap integral 

A" = gM 5 J dy (y; yf ) <p) (y; y]) <p h (y; y h ) , (4) 

with \^,v, M determined analogously. 2 In this letter the 
coupling g is assumed to be universal, as could result 



from higher-dimensional gauge interactions. We assume 
very small neutrino masses are due to the see-saw mech- 
anism. After integrating out the right-handed neutrinos 
V, this generates a low-energy left-handed Majorana neu- 
trino mass matrix from the effective interaction 

(C,y/ (</>)) k Ijhh, where Cy = {\ T U A M * \ v ).. . (5) 

To form the landscape we assume that the center coordi- 
nates yf of all wavefunctions arc scanned freely and inde- 
pendently of one another on S 1 . All of the other param- 
eters are set by hand. 3 This random scanning of center 
coordinates is motivated by some knowledge of instanton 
moduli; the center coordinates of multi-instantons can be 
chosen freely, and zero modes tend to localize around the 
centers of these instantons. 4 

Fig. [T] shows distributions of observables that follow 
from the Gaussian landscape on S . Yukawa eigenval- 
ues arc hierarchical (i) and all three mixing angles in the 
quark sector are small (ii), if the wavefunctions of the 
Higgs boson and fields in the 10 = (q, u, e) multiplets are 
localized in the extra dimensions, i.e., dio, dn <S L. The 
overlap between a quark doublet qi and the Higgs boson 
introduces a correlation between the up-type and down- 
type Yukawa matrices, giving the observed electroweak 
pairing (i.e. no peaks at 7r/2 for the mixing angle dis- 
tributions) in the quark sector (iii). Thus, the localized 
wavefunction of the Higgs boson, as well as those of quark 
doublets, are the essence of what wc observe as the gen- 
eration structure in the quark sector. We also find that 
the observed hierarchy 6»f 3 KM < 6»f 2 KM , 0§, KM is typical. 

The leptonic mixing angles arc typically of order unity 
(iv), when the width parameter d$ is set comparable to 
the size of the extra dimension, i.e. when the fields in 
5 = (d, I) do not have very localized wavefunctions. The 
angle #13 tends to be smaller than the other two neu- 
trino oscillation angles, agreeing very well with observa- 
tion (v). The larger d% width implies a milder hierar- 
chy among the Yukawa eigenvalues in the charged lep- 
ton and down-quark sectors than in the up-quark sector 
(vi). Although the distribution of CP phases depends on 
the details of how complex phases are introduced in the 
landscape, for the wavefunction ([2]) with r = 3 there are 
continuous and flat components in the distributions of 
Sckm and 6 V (vii). We consider it unlikely that these flat 



1 These are dual to the complex structure moduli of the Type IIB 
string theory (F-theory) in cases with J\f = 1 supersymmetry. 

2 Note that we use ip h (y;y h ) also for the overlap integrations for 
A^j e , which allows for the analytical understanding described 
later. Since the wavefunction of the Higgs boson is not scalar 
valued if it originates from a vector field in extra dimensions, the 
"wavefunction" ip h that is expected to be used for X d ' e does 
not have to be the complex conjugate of ip h . We emphasize that 
using a universal value for r is not more than one of the simplest 
ways to introduce complex phases into Gaussian landscapes. 



3 By slicing a subset out of a possibly much larger landscape of 
vacua, we see that there exists a subset that is phcnomcnologi- 
cally successful. It is a separate question whether such a subset is 
highly weighted in the vacuum statistics of the entire string land- 
scape or because of cosmological evolution and / or environmental 
selection. Wc consider that practical progress can be made by 
splitting the full problem into simpler parts. 

4 This is certainly not a rigorous argument; indeed basis inde- 
pendence of the probability distributions is also lost here (see 
Ref. Q). It would be interesting if this assumption — the ran- 
dom scanning of center coordinates — were refined by studying 
flux compactification. 
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FIG. 1: Distributions of observables in the S 1 Gaussian landscape. The width parameters are set at dh = dio = d^ — = 0.08L 
and = 0.3L, and we use r — 3 and g = 0.2. Numbers in brackets are the experimentally measured values (or limits), with 
the leading renormalization effects up to the Planck scale partially taken into account. All logarithms are base ten. 



components will disappear for minor modifications to the 
form of the wavefunctions. 



Neutrino masses generated via the see-saw mechanism 
are typically very hierarchical in Gaussian landscapes, 
arising from both the overlap between 77 and h in A 17 , 
and the overlap between V and <ft in \ M . These two hi- 
erarchies are not correlated unless the Higgs boson and 
(j) wavefunctions are correlated. Therefore, in the statis- 
tical distribution of see-saw neutrino mass eigenvalues, 
these hierarchies add. An inverted hierarchy is quite un- 
likely. We note that complex phases play a crucial role 
in allowing for a sufficiently large value of / m% Q • 

The S 1 Gaussian landscape does not very well explain 
the large observed values of At/Aft >T . Ref. Q discusses 
how this problem can be addressed in Gaussian land- 
scapes on geometries other than S 1 . 



3. Analytic Approximation: Most of the distribution 
functions in Fig. [T] can be understood with an analytic 
approximation. Let us begin with the down-type Yukawa 
matrix. The approximation is to consider d 10 ,dh <C L, 
with the various wavefunctions peaked not too far from 



each other, and as O(L). Then 



A % oc g <pf 



V = (dly q 1 +dl y h )/(dl 



ef)=exp [-(l + ri)(y]-y»f/2(dl 



4); vt 

\- dl)} . 



(6) 



The wavefunctions of the dj's (C 5) are not localized, yet 
not absolutely flat over S 1 . With the random scanning of 
the center coordinates yf , y\ and y h , the first factor <j>f ef- 
fectively yields random coefficients of order unity. On the 
other hand, the second factor can be exponentially 

small. The quantity In I e'W I = -(Ay) 2 /[2(d 2 10 +d%)] can 
be as small as — £ 2 /[8(c^ -r-<i^)] in the S 1 Gaussian land- 
scape, and its distribution function is determined only 
from the geometry of S 1 : dP/d\Ay\ = 2/L, such that 



dP 



'2(d? t 



dl) 



d\ln\eiW\\ I/ £ 2 |ln|e9( d )|| 



/(|ln|e^)||). (7) 



When a down-type Yukawa matrix is generated in the 5* 1 



Gaussian landscape, the three ln|e| |'s randomly follow 
a distribution function ([7]), and the largest, middle and 
smallest among them determine the order of magnitude 
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of ln(A{,/<7), \n(X s /g) and ln(Ad/<?), respectively. Thus 



d 3 P 



d\hi At | d|ln A s | e?|ln A^ 



3!/(|lnA & |)/(|lnA s |)/(|lnA d | 



.(8) 

for .g = l. The distribution of each of In \b, s ,d is obtained 
by integrating the two other variables. The resulting an- 
alytic approximate distribution agrees very well with the 
results of the numerical simulation above. 

The up-type Yukawa matrix has the structure 



v 



9ijt 



<?(«) 



(9) 



for the widths d q = d a = dio, with statistically neutral 
random coefficients g^j and "flavor suppression factors" 
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The e"'s are given by the same expression except that 
y q - is replaced by yf. These flavor suppression factors 
follow a distribution function similar to ([7]), and the dis- 
tribution functions of the largest, middle and smallest 
suppression factors are obtained just as in The dis- 
tribution function of In At (lnA c and lnA u ) is obtained 
by convoluting the distribution functions of the largest 
(middle and smallest) lne" and lne 9 ^"'. 

Similarly approximate distribution functions can also 
be derived for the mixing angles of the quark sector and 
the mass eigenvalues of the charged leptons and see-saw 
neutrinos. From the the flavor structure of the Yukawa 
matrices © and hierarchical quark masses and small 
mixing angles in the quark sector follow from Gaussian 
landscapes, just as in the conventional flavor symmetry 
approach. However, in Gaussian landscapes the geom- 
etry of the extra dimensions {dP/d{Ay) 2 in particular) 
automatically determines the "flavor symmetry charge" 
assignments within a given representation. Note also that 
the flavor suppression factor associated with q.j in 

the up-type Yukawa matrix is not the same as e? in 
the down-type Yukawa matrix. 

4. Geometry Dependence: So far we have consid- 
ered only the Gaussian landscape on S . Ultimately we 
would like to understand Gaussian landscapes with more 
extra dimensions and on non-trivial geometries. For ex- 
ample, three-dimensional base manifolds B, along with 
a T 3 fiber, can be used as six-dimensional manifolds in 
string theory compactification. Therefore B of more than 
one dimension will be of practical interest. Such a study 
is launched in Ref. 0] ; here we summarize initial results. 

The flavor structures (i)-(vi) observed in nature are 
obtained statistically in Gaussian landscapes for various 
geometries B, as long as Gaussian wavefunctions on B 
are localized for 10's and the Higgs but not for 5's. The 
distribution functions of obscrvables certainly depend on 
the geometry of B, but in the analytic approximation, 
they depend only on the volume distribution dP/d{Ay) 2 



of the geometry. The distribution functions of the obscrv- 
ables are obtained by multiplying and integrating the vol- 
ume distribution functions many times. Because of these 
integrations, the details in the volume distribution func- 
tion dP/d{Ay) 2 are smeared and only averaged proper- 
ties of the geometry of B are reflected in the distribution 
of observables. Some observables are more sensitive to 
the geometry, some not. See Ref. Q for more details. 

5. Conditional Probability: We have so far discussed 
vacuum statistics, assigning equal weight to each vacuum 
in the landscape. Questions of real interest, however, 
are probability distributions of obscrvables with various 
weight factors from cosmological evolution and environ- 
mental selection included. Furthermore, we are inter- 
ested in the distributions of yet-to-be measured observ- 
ables in the subset of vacua that pass all known experi- 
mental constraints. We hope to measure one mixing an- 
gle and one CP phase 8 U in future neutrino oscillation 
experiments, and one mass parameter mpp in neutrino- 
less double beta decay experiments. The leading contri- 
bution to mpp is U 2 2 rri2- mpp = U 2 2 m%(l + AR), but a 
precise prediction can be made only if the distribution for 
|Ai?| = \{U e i/U e 2) 2 {m\/m2)\ is peaked at small values. 

It is not practical to obtain precise estimates of the 
weight factors or to generate statistics large enough so 
that all of the experimental cuts can be imposed. In- 
stead, we use the Gaussian landscape on S 1 and ex- 
tract a subset that satisfy "loose experimental cuts" A) 
10" 2 < Am%/Am 2 atm < 1CT 1 , B) sin 2 (26» Q ) > 0.7, C) 
sin 2 (26» atm ) > 0.8 and, D) sin6»i 3 < 0.18, to get a feel- 
ing for how much the distributions of #13, 8 U and |Ai?| 
are affected by the weights and/or cuts. The results are 
displayed in Fig. [TJ We see that imposing these cuts 
makes a significant difference in the expectations for fu- 
ture experiments. In particular the distribution for #13 
becomes peaked near the experimental limit, and there 
is no longer a contribution to the 5 V distribution that is 
peaked around zero. Finally, the distribution for AR is 
greatly reduced for all values of AR above 0.1, sharpen- 
ing the prediction for mpp on the S 1 landscape. 

6. Conclusion: We have shown that many features of 
flavor can be broadly understood from a simple Gaus- 
sian landscape using just five free parameters. Distribu- 
tions for the 17 measured flavor parameters, as well as 
for neutrino parameters that may be measured in future 
experiments, are shown in Figure 1. The distributions 
are broad and result from scanning the centers of each 
quark, lepton and Higgs wavefunction randomly and in- 
dependently. Even if flavor originates from a landscape of 
vacua, it may be that some flavor observables are affected 
by environmental selection, for example the first gener- 
ation masses m u ^d,e may be selected for certain atomic 
and nuclear properties. 
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